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Abstract 

The non-local thermodynaic equilibrium (non-LTE) line formation for the two ions of 
zirconium is considered through a range of spectral types when the Zr abundance varies 
from the solar value down to [Zr/H] = —3. The model atom was built using 148 energy 
levels of Zr I, 772 levels of Zr II, and the ground state of Zr III. It was shown that the main 
non-LTE mechnism for the minority species Zr I is ultraviolet overionization. Non-LTE 
leads to systematically depleted total absorption in the Zr I lines and positive abundance 
corrections, reaching to 0.33 dex for the solar metallicity models. The excited levels of 
Zr II are overpopulated relative to their thermodynamic equilibrium (TE) populations in 
the line formation layers due to radiative pumping from the low-excitation levels. As a 
result, the line source function exceeds the Planck function leading to weakening the Zr II 
lines and positive non-LTE abundance corrections. Such corrections grow towards lower 
metallicity and lower surface gravity and reach to 0.34 dex for log g = 2.0 and [M/H] 
= —2 at T e fj = 5500 K. As a test and first application of the Zr I-Zr II model atom, Zr 
abundance was determined for the Sun on the basis of ID LTE model atmosphere. Lines 
of Zr I and Zr II give consistent within the error bars non-LTE abundances, while the 
difference in LTE abundances amounts to 0.28 dex. The solar abundance of zirconium 
obtained with the MAFAGS solar model atmosphere is logezr,© = 2.63±0.07. 



1 Introduction 

One believes that the elements beyound the iron group are produced by neutron-capture reac- 
tions. These reactions are subdivided into rapid (r) and slow (s) processes depending on the 
neutron flux available. The main and secondary or weak component are distinguished among the 
slow processes. The r-process is often associated with type-II supernova explosions. The weak 
and main components of the s-process can take place, respectively, in the core of high-mass 
(M > 20M & ) stars at the hydrostatic helium core burning stage and during the thermally- 
pulsing AGB phase of intermediate-mass (2 — 4M Q ) stars. According to Kappeler et al. [16J, 
in high-mass stars the nuclei of only the lightest of the heavy elements, with atomic masses 
A < 90, can be produced in the s-process. Since the atomic mass of zirconium is A = 91.22, it 
can be produced in all three types of neutron-capture reactions. The characteristic production 
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time is different for different types of reactions. Therefore, the relative contributions from each 
of the processes to the abundance of a specific element changed as the Galaxy evolved. At 
present, the theories of the r-process and the weak component of the s-process cannot accu- 
rately predict the yields of elements. Therefore, it is very important to restore the history of 
heavy elements enrichment of the interstellar medium and, thus, to impose constraints on the 
nucleosynthesis theories on the basis of observational data. 

In our previous paper |24| . we determined the zirconium abundance assuming the local 
thermodynamic equilibrium (LTE) for a sample of 52 stars belonging to three types of Galactic 
population - the thin disk, the thick disk, and the halo. We found a large overabundance of 
zirconium relative to barium, up to log(Zr/Ba) = 1.9 at [Ba/H] = —3.8, in the halo stars. 
The Zr/Ba ratio decreases towards higher barium abundance and reaches its solar value of 
log(Zr /Ba ) = 0.41 j20] in the thin-disk stars. If the heavy elements would synthesized only 
in the r-process in the early Galaxy, until the first intermediate-mass stars evolved, the ratio 
log(Zr/Ba) in old stars must be constant at a level of 0.33 dex [3J. Our observational data 
appeared to infer a distinct production mechanism for the light trans-iron (Sr-Zr) and heavy 
elements beyond Ba in the early Galaxy. Similar conclusions were drawn in [2J [T2] . 

In this paper, we want to check whether departures from LTE can affect the determined 
evolutionary behavior of the zirconium abundance. We develop a technique for calculating the 
statistical equilibrium (SE) of Zr I-Zr II and consider the non-local thermodynamic equilibrium 
(non-LTE) line formation for Zr I and Zr II in the atmospheres of the Sun and cool stars. 

The paper is structured as follows. In the next section, we describe the Zr I-Zr II model 
atom, the method of non-LTE calculations, and the mechanisms of departure from LTE for 
Zr I and Zr II. In Sect.[3l we analyze the solar zirconium lines, determine the solar zirconium 
abundance from Zr I and Zr II lines for various line formation scenarios, and compare our LTE 
results with the literature data. The non-LTE effects for Zr I and Zr II lines depending on 
stellar parameters are presented in Sect.HJ Finally, we give our conclusions. 



2 Non-LTE calculations 
2.1 Model atom 

Energy levels. In the atmosphere of cool stars with effective temperature T e g > 5000 K, zir- 
conium is largely ionized and the fraction of Zr I atoms does not exceed a thousandth of the 
total number of zirconium atoms. Therefore, the Zr I level populations are very sensitive to 
a change in the net ionization minus recombination rate. The Zr II ionization stage is of our 
particular interest, because only Zr II lines are detected in the metal-deficient stars. 

The system of Zr I measured levels includes the singlet, triplet, and quintet terms of the 
electronic configurations 4d 2 5s 2 ; 4d 2 5s5p; 4d 3 nl, where nl = 5s and 5p, in total, 148 levels with 
excitation energy up to E exc = 4.25 eV from the NIST (http : / / physics .nist. gov/1 PhysRefData) 
and VALD [17] atomic data databases. The ground state of Zr I is 4d 2 5s 2 3 F2. 

The experimental Zr II levels (a total of 103) provided by the NIST database and by [23J 
have E exc < 8.75 eV and belong to the doublet and quartet terms of the electronic configurations 
4d 2 nl (nl = 5s, 5p), 4d5s 2 , and 4d 3 . The ground state of Zr II is 4d 2 5s 4 F 3 / 2 . We also employ 
669 levels predicted in our calculations of the Zr II atomic structure. These levels belong to 
the doublet and quartet terms of the following electronic configurations: 4d 3 , 4d 2 nl, where nl 
= 5p, 5d, 5f, 5g, 6s, 6d, 6p, 7s, 7d, and 7p, and 4d5snl, where nl = 4f, 5p, 5d, 6s, 6p, 6d, and 
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Figure 1: The term diagram of Zr I-Zr II as known from laboratory measurements and atomic 
structure calculations for Zr II. 

7p. The term diagram of Zr I and Zr II is shown in FigJTJ 

In constructing the model atom, the fine structure was taken into account for all the Zr II 
levels with E exc < 5 eV. The remaining levels with common parity and close energies were 
combined. The final model atom includes 63 levels of Zr I, 247 levels of Zr II, and the ground 
state of Zr III. 

Radiative bound-bound (b-b) rates. The system of the Zr I levels in our model atom includes 
246 permitted bound-bound transitions. The oscillator strengths (fij) were taken from the 
NIST and VALD databases. The Zr II model atom includes 9336 radiative b — b transitions. 
For 1070 of them, we use the experimental values from [22j and the data from the VALD 
database. For the remaining transitions, we apply /jj-values calculated in this study. 

Photoionization cross sections. Since there are no accurate data on the photoionization cross 
sections for Zr I and Zr II, we apply a hydrogenic approximation with an effective principal 
quantum number n e ff\ 

n e ff = Z eff 

where Z e ff is the effective atomic charge (Z e // = 1 for neutral atoms, Z e ff = 2 for the first 
ions, etc.), \h is the ionization energy of the hydrogen atom from its ground state, and \ is 
the level ionization energy. 

Collisional rates. The electron impact excitation was taken into account using the formula 
from |26| for permitted transitions. For forbidden transitions, the effective collision strength 




was assumed to be T = 1. The electron impact ionization was taken into account using the 
formula from In the atmosphere of cool stars, the number density of neutral hydrogen 

atoms is much greater than that for electrons. Therefore, it is also important to take into 
account the inelastic collisions with neutral hydrogen atoms. We used the Steenbock & Holweger 
formula from [27J for permitted transitions and the approximation of Takeda [2B] for forbidden 
transitions. Since both formulas provide an accuracy only in order of magnitude, we introduced 
a scaling factor in the range from to 1, which we attempted to determine by analyzing the 
solar Zr I and Zr II lines. 



2.2 Mechanisms of departure from LTE 

In all our calculations, we used plane-parallel, homogeneous, blanketed model atmospheres 
computed using the MAFAGS code [13]. The combined system of radiation transfer and statis- 
tical equilibrium equations was solved using the DETAIL code [10] based on the accelerated A 
iteration method. The level populations obtained were used to calculate the profiles of the lines 
under study. Synthetic spectra were computed using the SIU code developed at the Munich 
University [25]. The list of spectral lines includes all atomic and molecular lines from the tables 
by Kurucz [19J. 

The departures from LTE in the level population are characterized by b— factor: 6j = rij/nj*, 
where n« is the level population obtained from the SE equations (non-LTE population) and rii* 
is the LTE population calculated from the Boltzmann-Saha formulas. Figure [2] shows the 
b— factors of some Zr I and Zr II levels in the solar atmosphere. The ground state and low 
excitation levels of Zr II keep their TE populations almost over the entire atmosphere. Only 
in the uppermost layers with logrsooo < ~ 3, the low excitation levels are slightly underpopu- 
lated, while the ground state is overpopulated. All of the remaining levels are overpopulated 
throughout the atmosphere, except for the 4d 2 5p 4 F^ 2 level with E exc = 3.69 eV at depths 
around logrsooo — — 3. The overpopulation of high-excitation Zr II levels results from radiative 
pumping of the ultra-violet (UV) transitions arising from the ground state and low-excitation 
levels. The main non-LTE mechanism for Zr I is the overionization caused by superthermal 
UV radiation of a non-local origin below the thresholds of the low excitation levels. All Zr I 
levels are underpopulated upwards logr 50 oo ~ 0.5. 

The Zr II lines used in the abundance determinations are formed in the solar atmosphere 
between logrsooo = and logrsooo = — 2. In these layers, the upper level of each transition is 
overpopulated relative to its LTE population, while the populations of the lower levels keep 
their TE values. The line source function is defined as 



13 ~ c 2 h fex V {hvlkT)-\' 

where bj and 6j are the b— factors of the upper and lower levels, respectively. For the lines in 
the visual spectral range (hu > kT), the following approximate expression is valid: 

where B V (T) is the Planck function. In the line formation layers, bj > bi is valid for each Zr II 
line (Fig. [2]), and the line source function is greater than the Planck function. As a result, 
the line is weakened compared to LTE and the zirconium abundance in non-LTE calculations 
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Figure 2: Departure coefficients log b for selected levels of Zr I and Zr II as a function of optical 
depth logT5ooo i n the solar model atmosphere. Departure coefficients of the ground states of 
Zr I and Zr II are plotted by continuous curves. 



should be increased to describe its LTE profile and equivalent width. Thus, the non-LTE effects 
lead to positive non-LTE corrections for Zr II lines. 

Non-LTE leads to a weakening of the Zr I lines, too. However, the mechanism of departures 
from LTE is different from that for Zr II lines. The Zr I lines are weaker relative to their 
LTE strengths mainly due to the general overionization and also due to bj/b{ > 1 resulting in 
Siu > B v and the depleted line absorption. At the depths where the Zr I lines are formed, the 
upper levels are all depleted to a lesser extent relative to their LTE populations than are the 
lower levels. The non-LTE abundance corrections are positive for the Zr I lines. 



3 Analysis of the solar zirconium lines 
3.1 Selection of lines 

Ljung et al. [22] and Biemont and Grevesse [9] provided the lists of lines which were used to 
determine the zirconium abundance of the solar atmosphere. We checked all lines from these 
lists for the presence of blends whose influence was difficult or impossible to take into account. 



For this purpose, we computed synthetic spectra using the SIU code by taking into account all 
of the atomic and molecular lines from the list by Kurucz [19] and compared them with the 
observed solar spectrum in fluxes [18] with a resolution R = 340 000 in the range 4000-4700 A 
and R = 500 000 at longer wavelengths. We use a solar model atmosphere with T e ff = 5780 K, 
logg = 4.44, and [M/H] = 0.0. The microturbulence velocity was assumed to be fixed, V m j C = 
0.9 km s _1 . To describe the observed profiles, the theoretical line profiles were convolved with a 
profile that combines a rotational broadening of 1.8 km s -1 and broadening by macroturbulence 
with a radial-tangential profile of 3 km s _1 to 4kms~ 1 for different lines. 

Tables [TJ and [5] list all of the checked Zr I and Zr II lines. The oscillator strengths of Zr II 
lines were taken from [22j as the most recent ones and from [9]. The data for Zr I lines are 
available only in j9]. All were obtained from laboratory measurements. Below, we give a 
brief description of all lines. 

Zr I lines. Biemont and Grevesse [9] used 34 Zr I lines to determine the zirconium abun- 
dance. It emerged that almost all of the lines from this list, except two lines, Zr I 4241.706 and 
4687.799 A, were strongly blended and we excluded them from the abundance determinations. 
Let us illustrate this with several examples. The strongest Zr I 3891.383 A line lies in the wing 
of the 3889.051 A line and is described without confidence. The Zr I lines are mostly very 
weak - only four of them have an equivalent width larger than 10 mA. Therefore, they can 
be affected even by weak lines of other atoms and molecules. Figure [3] shows our attempts to 
describe the blends near four Zr I lines. Zr I 4710.077 A is one of the strongest lines, but it 
is located in a molecular band. Besides, it is influenced by two strong lines, Ti I 4710.183 A 
and Ti I 4710.189 A. The Zr I 3509.331 A line is influenced by the CN 3509.247 A and Ni l 
3509.362 A lines and by the absorption of an unknown origin in the range 3509.2-3509.6 A. 
We are unable to properly take into account the absorption around the Zr I 3601.198 A line, 
because we do not know what is responsible for it. The Zr I 4030.049 A line is weak (W\ 
= 2.6 mA, Table [TJ and is lost in the absorption produced by the wings of the strong Fe I 
4030.185 A line. 

The profiles of the strong Zr I 4043, 4072, 4739, 4772 A lines are also distorted by the 
influence of molecular and atomic lines. In addition, in most cases, there is absorption near 
the zirconium lines that we are unable to take into account, because we do not know what is 
responsible for it. The profiles of the lines being studied are generally asymmetric, suggesting 
the presence of other absorption sources at close wavelengths. 

Zr II lines. We checked 25 Zr II lines, 16 of which turned out to be unsuitable for the 
abundance determination. Let us illustrate this with examples. Figure H] shows examples of 
the Zr II lines excluded from the subsequent analysis. The Zr II 3458.932 A is located in the 
wing of the strong Ni I 3458.456 A line and is also blended with Ti II 3458.903 A. As can be 
seen from Fig. HI we cannot describe well the local continuum near the Zr II 3458.932 A line. 
The Zr II 3454.572 A line is blended with Fe I 3454.569 A and we cannot reproduce the local 
continuum. We cannot describe the absorption near Zr II 3432.404 A either, because we do not 
know was is responsible for it. 

The Zr II 4024.435 A line from the list by Ljung et al. [22J deserves particular attention. 
Despite good agreement with the abundances determined by Ljung et al. j22] from different 
lines, we believe that it cannot be used in our analysis because of strong blending. The Zr II 
4024.435 A line is influenced by three lines at once: Ce II 4024.485 A, Fe II 4024.547 A, and 
Ti I 4024.572 A (Fig. H]). To achieve the best agreement between the observed and theoretical 
spectra, the oscillator strengths of the cerium and titanium lines were reduced by 0.07 and 
0.10 dex, respectively. The figure shows the contribution from each of the lines to the total 



Table 1: Atomic data for the Zr I lines, solar equivalent widths, and LTE abundances derived 
in this study and by Biemont & Grevesse [9]. 



A s , A Xexc, eV \oggf B Wf , mA \oge z r transition/comment 



4241.706 0.65 0.14 3.7 2.33 y 4d 3 5s 5 F 3 - 4d 3 5p 5 F° 

4687.805 0.73 0.55 10.0 2.34 y 4d 3 5s 5 F 5 - 4d 3 5p 5 G° 6 
Lines excluded from zirconium abundance determination 
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Figure 3: Solar Zr I lines excluded from the abundance determination. Synthetic LTE (continu- 
ous curve) flux profiles of the 3509, 3601, 4030, and 4710 A blends compared with the observed 
spectrum of the Kurucz et al. [18] solar flux atlas (bold dots). Dash-dotted curve shows the 
pure Zr I 4710 A line, while dotted curve shows the molecular C2 lines. Dashed curve plots the 
synthetic spectrum of the 3509 A blend calculated with no zirconium in the atmosphere. 

absorption. We exclude the Zr II 4024.435 A line from the subsequent analysis. 

The Zr II 3549, 3588 A lines are located in molecular bands. The Zr II 3479.017, 3499, 
3607, 3671, 3714, 3796, 3836, 4034, 4085, 4317 A lines are strongly influenced by the lines of 
other atoms and the accuracy of determining the zirconium abundance depends strongly on 
the accuracy of the atomic parameters of the blending lines. Below, we give the list of Zr I and 
Zr II lines used to determine the zirconium abundance with the corresponding comments. 

Zr I 4241-706 A. As can be seen from Fig. [51 the Zr I 4241.706 A line is influenced by the 
SiH molecular line (W\ = 0.6 mA, 11 of the total equivalent width). As /y(SiH) increases by 
0.2 dex, the abundance determined from the Zr I 4241.706 A line decreases by 0.02 dex. The 
total profile is described more poorly. 

Zr I 4687.805 A. This line is blended with the Ti I 4687.809 A line (W A = 1.1 mA) and 
the right wing of the Zr I 4687.805 A line being studied is influenced by the C2 4687.9 A 
molecular line, as shown in Fig. |5j Changing /y(C2) by 0.5 dex does not lead to any change in 
the zirconium abundance determined from the Zr I 4687.805 A line. The atomic data for Ti I 
4687.809 A were taken from VALD. 

Zr II 3479.387 A. As can be seen from Fig. El the observed wings of this line lie below 
the theoretical ones approximately by 5%, suggesting that we do not completely take into 
account the absorption in the local continuum. When the local continuum was raised by 5%, 
the equivalent width of the Zr II 3479.387 A line decreased by 6.4%. In this case, the zirconium 
abundance decreased by 0.13 dex. 

Zr II 3505.666 A. The theoretical level of the local continuum near this line is above the 




Figure 4: Solar lines of Zr II excluded from the abundance determination. Synthetic LTE 
(continuous curve) flux profiles of the 3432, 3454, 3458, and 4024 A blends compared with the 
observed spectrum of the Kurucz et al. [18j solar flux atlas (bold dots). Dash-dotted curve 
shows the pure Zr II 4024 A line, while long-dashed curve shows the lines of Ce II, Fe II, 
and Ti I. Dashed curve plots the synthetic spectrum of the 3454 A blend calculated with no 
zirconium in the atmosphere. 
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Figure 5: Synthetic LTE (continuous curve) flux profiles of the selected Zr I lines compared 
with the observed spectrum of the Kurucz et al. [18] solar flux atlas (bold dots). Dash-dotted 
curve shows the pure Zr I 4241 and 4687 A lines. Dotted curve shows the molecular SiH and 
C2 lines. Long-dashed curve shows the Ti I line. 



observed one by 2%. Since the line is strong (Wa = 53.4 mA), this affects insignificantly the 
abundance determination from this line. When the local continuum is raised by 2%, the line 
equivalent width W\{7^i II) decreases by 3.7% and the abundance decreases by 0.07 dex. 

Zr II 3551.951 A. This line is located in the CN molecular band whose influence is weak 

(Fig. ED- 

Zr II 4050.320 A. This line is blended with two lines of manganese, Mn I 4050.288 and 
4050.345 A (Fig. [6]), whose abundance in the solar atmosphere is known. Therefore, we use 
the Zr II 4050.320 A line in analyzing the zirconium abundance only in the solar atmosphere. 
However, it should be used with caution for other stars and only when the manganese abundance 
is determined carefully. 

Zr II 4208.980 A. The red wing of this line is slightly influenced by the CN 4209.046 A 
molecular line. According to the list by Kurucz p9], the equivalent width of the molecular line 
should not exceed 1.4 mA, while W\(Zy II 4208 A) = 44.4 mA. Figure [6] describes the observed 
spectrum near the Zr II 4208.980 A line with allowance made for all lines and separately shows 
the contribution from the molecular line. Increasing logp/(CN) by 0.1 dex causes the zirconium 
abundance determined from the Zr II 4208.980 A line to decrease by less than 0.01 dex. 

Zr II 4^58.041 A. The stronger Fe II 4258.154 A line affects only the red wing of the 
zirconium line (Fig. [6]). To describe well the iron line, we had to increase its wavelength by 
0.06 A and the oscillator strength by 0.28 dex. However, this increase in logg/(Fe II) had 
virtually no effect on the abundance determined from the Zr II 4258.041 A line (A(log£z r ) < 
0.01 dex). 

Zr II 444^-992 A. The zirconium line lies between two strong lines, Fe I 4442.832 and 
4443.194 A (Fig. [6]). For the best description of the observed spectrum, we had to increase the 
oscillator strength of the Fe I 4442.832 A line by 0.05 dex. Increasing the oscillator strengths 
of both iron lines by 0.05 dex causes the zirconium abundance determined from the Zr II 
4442.992 A line to decrease by 0.01 dex. 

Zr II 4496.962 A. This line is located in the wing of the Cr I 4496.854 A line (Fig. EJ. 
We determined the zirconium abundance from the red wing of the Zr II 4496.962 A line. 
The oscillator strength had to be increased by 0.1 dex, which had no effect on the zirconium 
abundance determination. 

Zr II 5112.210 A. This line is virtually isolated. The C2 5112.361 A molecular line makes a 
small contribution to the absorption. Changing logp/(C2) by 0.1 does not lead to any change 
in the zirconium abundance determined from the Zr II 5112.270 A line. 

We concluded that most of the lines from the list by Biemont and Grevesse [9] and one Zr II 
line from the list by Ljung et al. [22J cannot be used to determine the zirconium abundance, 
because the contribution from neighboring lines cannot be properly taken into account, espe- 
cially when the absorption near the line under study is produced by unknown sources (i.e., the 
lines are absent in the list by Kurucz [19J). Only two Zr I lines and nine Zr II line can be used 
to reliably estimate the zirconium abundance in the solar atmosphere. 

Zirconium in nature is represented by five isotopes: 90 Zr, 91 Zr, 92 Zr, 94 Zr, and 96 Zr. In 
the matter of the Solar system, the relative contribution from each of the isotopes to the total 
zirconium abundance is 51.45:11.22:17.15:17.38:2.80 [20J. None of the works on the experimental 
study of the of the zirconium spectrum provides data on the isotopic shifts and hyperfine level 
splitting. Therefore, we considered the zirconium lines to be one-component ones. 



Table 2: Atomic data for the Zr II lines, solar equivalent widths, and LTE abundances derived 
in this study, Biemont & Grevesse [9j, and Ljung et al. [22J. 
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Figure 6: Synthetic LTE (continuous curve) flux profiles of the selected Zr II lines compared 
with the observed spectrum of the Kurucz et al. [18J solar flux atlas (bold dots). Dash-dotted 
curve shows the pure Zr II 4258, 4443, and 4496 A lines. Dotted curve shows the molecular 
CN lines. Long-dashed curve shows the Mn I lines. 



3.2 Solar zirconium abundance 



Using the selected Zr I and Zr II lines, we determined the mean zirconium LTE abundances for 
each ionization stage (Table [3]). The difference loge(Zr II) - loge(Zr I) = 0.28 dex turned out 
to exceed considerably the error in the abundance determined from individual lines. 



Table 3: The solar Zr abundances from the individual Zr I and Zr II lines and the mean values 
derived from the LTE and NLTE (k H = 0, 0.1, 0.33, 1, 0.1+, 0.33+, and 1+) calculations with 
the MAFAGS solar model atmosphere of T cff = 5780 K, hgg = AAA, V mic = 0.9 km/s. 



version 



Zr I 
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Zr II 
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NLTE 

k H =0.00 

k H =0.10 

k H =0.33 

k H =l.QQ 

k H =0.l0+ 

k H =0.33+ 

fc H =1.00+ 



2.33 2.34 2.33 0.01 

2.68 2.61 2.65 0.05 

2.68 2.61 2.65 0.05 

2.67 2.59 2.63 0.06 

2.65 2.57 2.61 0.06 

2.62 2.61 2.62 0.01 

2.57 2.58 2.58 0.01 

2.53 2.54 2.54 0.01 



2.66 2.63 2.49 2.42 2.65 2.61 2.65 2.67 2.67 2.61 0.09 

2.77 2.67 2.54 2.50 2.76 2.65 2.73 2.75 2.82 2.69 0.11 

2.77 2.67 2.54 2.50 2.76 2.65 2.72 2.74 2.81 2.68 0.11 

2.76 2.66 2.53 2.49 2.75 2.64 2.72 2.73 2.80 2.68 0.11 

2.76 2.65 2.52 2.48 2.73 2.64 2.71 2.72 2.79 2.67 0.11 

2.74 2.65 2.52 2.45 2.70 2.62 2.70 2.70 2.72 2.64 0.10 

2.72 2.64 2.51 2.43 2.68 2.62 2.67 2.69 2.69 2.63 0.10 

2.70 2.63 2.50 2.42 2.66 2.61 2.66 2.68 2.68 2.62 0.09 



When the LTE assumption is abandoned, the accuracy of the derived zirconium abundance 
depends on the accuracy of the photoionization cross sections <j p h, transition probabilities, and 
the rates of inelastic collisions with electrons, Cjj(e), and neutral hydrogen atoms, Cij(H). 
The influence of an uncertainty in a p h and Cjj(e) on the results will be considered in the next 
subsection. Here, we will talk about the collisions with hydrogen atoms and about the oscillator 
strengths. 

According to recent studies by Belyaev et al. [7] and Belyaev and Barklem |8] for the 
resonance Na I and Li I lines, the formula from Steenbock and Holweger [27J gives overestimated 
rates of collisions with hydrogen atoms. Therefore, we performed our calculations with a scale 
factor kn that was varied in the range from to 1. We made an attempt to estimate kn 
empirically by analyzing Zr I and Zr II lines in the solar spectrum. The calculations were made 
for eight spectral line formation scenarios: LTE, non-LTE without allowance for the collisions 
with hydrogen atoms, non-LTE with allowance for the collisions with hydrogen atoms only 
in permitted transitions with kjj = 0.1, 0.33, 1.0, non-LTE with allowance for the hydrogen 
collisions in both permitted and forbidden transitions with kn = 0.1, 0.33, 1.0 (below, we will 
denote them by 0.1+, 0.33+, and 1.0+). For each of the scenarios, we determined the mean 
abundances from the lines of two ionization stages, Zr I and Zr II. The results are presented in 
Fig. El 

The left and right panels of Fig. [7] differ only by the use of different sources of oscillator 
strengths for Zr II lines, i.e., |22j and [9j, respectively. On the left panel, there exists a systematic 
abundance difference between the two ionization stages: the zirconium abundances derived 
from Zr II lines are higher than those inferred from Zr I lines independent of the line formation 
assumptions used. It should be noted that the oscillator strengths of Zr II lines from Ljung 
et al. are systematically lower than those from and Grevesse [9J by a value from 0.03 dex 
for the Zr II 3505 A line to 0.26 dex for the Zr II 5112 A line. When the oscillator strengths 




Figure 7: The mean solar Zr abundances derived from the lines of Zr I (squares) and Zr II 
(triangles) using the LTE and NLTE (k H = 0, 0.1, 0.33, 1, 0.1+, 0.33+, and 1+) calculations. 
In left and right panel, we use oscillator strengths of the Zr II lines from Ljung et al. [22J and 
Biemont & Grevesse [9], respectively. The meteoritic Zr abundance from Anders & Grevesse 
PQ, Lodders et al. [20], and Asplund et al. |1] is indicated by dash-dotted, dotted, and dashed 
curves, respectively. 

from [9j are used, the mean abundances derived from Zr II lines decrease by 0.09 dex in all our 
non-LTE calculations and the systematic difference log£^ r (Zr II) — log£zr(Zr I) changes its 
sign (Fig. [7]) compared to the case where fij from Ljung et al.[22] are used for Zr II lines. This 
points to the need for revising the values of for Zr I lines. Below, we consider only the case 
where from Ljung et al.[22] are used for Zr II lines. We see that any of the non-LTE scenarios 
leads to better agreement from the lines of the two ionization stages than the LTE case: the 
abundances inferred from Zr I and Zr II agree between themselves and with the meteorite 
abundances within la. Since we cannot give preference to any of the non-LTE calculations, we 
abandon the attempt to empirically estimate Nevertheless, for the subsequent study, we 
dwelt on ku = 0.1+, for which the difference loge^ r (Zr II) — log£zr(Zr I) is minimal and equal 
to 0.02 dex (see Table |3]). 

3.3 Influence of the uncertainties in atomic parameters on the non- 
LTE results 

We performed test calculations by varying the photoionization cross sections a p h and the rates 
of collisions with electrons C^(e). We traced the changes in the abundance derived from all 
Zr I and Zr II lines when the photoionization cross sections increased and decreased by a 
factor of 10: a p h{n e ff) xlO and <j p / i (n e //)/10 (for the standard calculations with kn = 0.1+ and 
using n e ff in the formula for a p h). Similar calculations were made for kn = 0.1+ but using the 
principal quantum number n, i.e., with <J p h(n), <r ph (n)xl0, a ph (n)/10. The results are presented 
in Table H 

When the principal quantum number is used instead of n e ff, the non-LTE effects for Zr I 
lines are reduced: the correction to the abundance for both Zr I lines is 0.04 dex. This is 
because n > n e ff for all Zr I levels; therefore, a p h{n) < a p h{n e ff and the depopulation of Zr I 
levels is less intense when using n than that when using n e ff. Increasing a p h(n) by a factor 
of 10 causes the abundance determined from Zr I lines to increase by 0.01 dex. A great effect 



Table 4: Influence of an uncertainty in atomic parameters on the non-LTE zirconium abundance 
determination. 
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(-0.09 dex) is obtained when a p h is reduced by a factor of 10, but it leads to a reduction in 
the departures from LTE for Zr I and, hence, to poorer agreement between Zr I and Zr II. 
When using n e ff, increasing the photoionization cross sections u p h{n e ff) xlO does not lead to a 
noticeable change in the departures from LTE, while the abundance determined from Zr I lines 
with a p h{n e ff) /10 decreases by 0.01 dex. 

Varying the photoionization cross sections affects very weakly the Zr II level populations, 
because this is the dominant ionization stage. The change in the abundance determined from 
Zr II lines is everywhere smaller than 0.01 dex. We also calculated the change in abundance 
when varying the rate of collisions with electrons: CV,(e)xl0 for kfj = 0.0 and kfj = 0.1+ 
(Table H]). The influence of collisions with electrons on the statistical equilibrium of zirco- 
nium must be strongest in the absence of hydrogen collisions. With kfj = 0.0, the difference 
\ogezr{Cij{e) x 10 — Cij(e)) ranges between —0.01 and —0.05 dex for different Zr I and Zr II 
lines, except the Zr II 5112 A line for which this difference is —0.08 dex. However, when the 
hydrogen collisions are taken into account, the influence of an uncertainty in CV,(e) on the 
departure from LTE decreases, because Cij(H) is an order of magnitude larger than Cij(e) for 
the transitions between levels with a small energy difference at which the collisional processes 
are particularly important. For example, for the standard case of kn = 0.1+, the difference 
\ogezr{Cij{e) x 10 — Cij(e)) is from to -0.03 dex. The sign of the effect is the same for the 
lines of both ionization stages and, hence, the difference between the mean abundances changes 
by no more than 0.01 dex. 

We see that the errors in the zirconium abundance due to the uncertainty in the photoion- 
ization cross sections and in the rate of collisions with electrons are smaller than the error in 
the mean abundance attributable to the spread in values obtained from different lines (a = 
0.07 dex). 

Thus, allowance for the departures from LTE makes it possible to reconcile, within the error 
limits, the zirconium abundances determined from the lines of the two ionization stages and we 
obtain the mean zirconium abundance \ogEz r = 2.63+0.07 (for kn = 0.1+). 



Table 5: Equivalent width and abundance comparisons between this study and Ljung et al. [22J 
for the selected lines of Zr II. Everywhere, the calculations were performed with the Holweger 
& Muller [15] solar model atmosphere. 
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3.4 Comparison with other studies 

The LTE zirconium abundance on the Sun has been determined previously by several authors. 
One of the first studies is the work by Biemont and Grevesse [9J , where the authors used 34 Zr I 
lines and 24 Zr II lines and the method of equivalent widths. The mean abundances derived 
from the lines of each ionization stage agree well between themselves and with the meteorite 
abundance, but with large errors: log£^ r (Zr I) = 2.57 ±0.21 and loge^ r (Zr II) = 2.56 ±0.14. 
In our view, the large errors are attributable to the difficulty of allowance for the blending of 
most of the Zr I and Zr II lines used. 

We pay particular attention to the study by Ljung et al |22| . They determined the zirconium 
abundance from seven Zr II lines (Table [2]) by the synthetic spectrum method using the solar 
spectrum in intensities and a three-dimensional model solar atmosphere. Subsequently, the 
equivalent widths of Zr II lines were calculated from the theoretical spectrum and the zirconium 
abundances were also determined using the MARCS [H] and Holweger- Muller [15] models from 
equivalent widths. 

For comparison, we made additional calculations with the Holweger-Muller model and the 
solar spectrum in intensities from the disk center from the atlas by Brault and Testerman [6]. 
From the private communication by Asplund |5j, we know that Ljung et al. [22] used the 
microturbulent velocity V m j C = l.Okms -1 . Therefore, we also took V m j C = 1.0 km s -1 for our 
calculations. The results are presented in Table [5J The equivalent widths of Zr II lines that we 
measured in the intensity spectrum turned out to agree with those from Ljung et al.[22] within 
1 mA, except for the Zr II 4442.992 and 4050.320 A lines for which the difference between our 
values of WA and those determined by Ljung et al.|22j is 1.5 and 2.9 mA, respectively. In this 
case, our abundances are systematically higher than those obtained by Ljung et al. [22] by a 
value from 0.07 to 0.11 dex for all of the lines, except Zr II 4050.320 A for which the difference of 
our values and those obtained by Ljung et al. [22] is Ae^r = —0.04 dex. The mean abundance 
from our determinations, \ogez r = 2.68 ± 0.06, is higher by 0.05 dex that the value obtained 
by Ljung et al.|22], \oge Zr = 2.63±0.02. 

The difference between our values of WA and those of Ljung et al.[22] can, in part, be 
explained by the subtleties of line profile fitting. For example, according to the private com- 
munication by Asplund [5J, they corrected their equivalent width of the Zr II 5112 A line and 
obtained a new value of WA(5112) = 8.5 mA, which agrees with our determination. However, 



Table 6: The NLTE abundance corrections for the Zr I and Zr II lines depending on stellar 
parameters. 
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the zirconium abundance determined from this line changed by only 0.01 dex: the value from 
Ljung et al. [22] is \ogez r = 2.64, while our value obtained from this line is logezr = 2.76. Thus, 
the difference in line equivalent widths cannot explain the difference in derived abundances. 

The cause of the discrepancies in zirconium abundance between our determinations and 
those of Ljung et al. (22] remains unclear. It may be related to differences in the techniques of 
calculations. 

4 Non-LTE effects depending on stellar parameters 

Here, we investigate the influence of non-LTE effects on Zr I and Zr II lines for a small grid 
of model atmospheres with effective temperatures T e ff = 5500 and 6000 K, surface gravities 
logg = 2.0 and 4.0, and metallicities [M/H] = —3, —2, —1, 0. Everywhere, microturbulence 
velocity is V m ic = 1.0 km s -1 . In these calculations, we took into account the collisions with 
neutral hydrogen atoms in both permitted and forbidden transitions with kn = 0.1. The results 
are presented in Table [6] as differences between the non-LTE and LTE abundances: Anlte = 
log^NLTE - log^LTE- Hereafter, Anlte is referred to as non-LTE correction. The non-LTE 
corrections are positive in all cases for both Zr I and Zr II lines. 

The non-LTE effects are strong for the Zr I lines in the solar metallicity models, with Anlte 
reaching to 0.33 dex for T e ff = 5500 K and logg = 4.0. We do not give data for the lower 
metallicity models because of very small equivalent widths of the Zr I lines (W\ < 1 mA). An 
exception is Zr I 4687 A in the model 5500/4.0/ — 1. 

Figure [8] shows the non-LTE corrections for the Zr II 3479 and 3551 A lines depending 
on metallicity, effective temperature, and surface gravity. These lines can be detected in low- 




metallicity stars, down to [M/H] = —3. As can be seen from Fig. |S]and Table El the abundance 
corrections are small for Zr II lines at solar metallicity: Anlte < 0.08 dex. However, it should 
be remembered that disregarding the departures from LTE introduces a systematic error that 
underestimates the derived abundance. The departures from LTE for the Zr II 3551 A line 
increase with decreasing metallicity. For Zr II 3479 A in the models with logg = 2.0, the 
non-LTE corrections increase from 0.01 dex at [M/H] = up to ~0.34 dex at [M/H] = —2. 
This is explained by an increase in the ultraviolet flux, on the one hand, and by a decrease 
in the number density of free electrons, on the other hand, and, as a result, by an increase in 
the radiative rate and a decrease in the collisional rate. As the metallicity decreases further, 
the non-LTE correction drops to 0.2 dex at [M/H] = —3 due to shifting the line formation 
depth to deeper atmospheric layers where collisions are efficient to establish thermodynamic 
equilibrium. A similar behavior is also observed for the Zr II 3479 A line in the models with 
T e ff = 5500 K and logg = 4.0 and T e ff = 6000 K and logg = 4.0. The non-LTE correction 
increases from 0.06 dex and 0.08 dex at [M/H] = to 0.14 dex and 0.17 dex at [M/H] = -1 
and then decreases to 0.11 dex and 0.12 dex at [M/H] = —3. The departures from LTE for 
Zr II lines are weakly sensitive to changes in effective temperature, which can also be seen from 
Fig. El 



5 Conclusions 

We constructed a model zirconium atom that includes 63 levels of Zr I , 247 levels of Zr II, and 
the ground Zr III state. We showed that the Zr I levels are underpopulated relative to their LTE 
populations due to ultraviolet overionization, while the excited Zr II levels are overpopulated 
because of radiative pumping from the ground state and lower excited levels. Non-LTE leads 
to weakening both Zr I and Zr II lines relative to their LTE strengths and positive non-LTE 
abundance corrections. 

The departures from LTE are strong for the Zr I lines in the solar metallicity models, with 
Anlte at the level of 0.3 dex. For most Zr II lines in the models with logg = 4.0, Anlte < 
0.12 dex. An exception is Zr II 3479 and 5112 A, with non-LTE corrections of 0.17 and 0.16 dex, 
respectively. The departures from LTE for Zr II are predicted to be stronger in giant stars (logg 



= 2). For example, A NLTE (Zr II 3479 A) = 0.34 dex in the model with T e// = 5500 K, logg = 
2.0, and [M/H] = —2. All the stars studied by Mashonkina et al. [23] at the LTE assumption 
are dwarfs or subgiants, with logg > 3.12. Taking the non-LTE effects for Zr II into account 
would lead to an up to 0.2 dex higher Zr overabundance relative to Ba than that found in [24j. 
Thus, we support our earlier conclusion about the growth of the Zr/Ba ratio towards lower 
barium abundance. 

From our analysis of Zr I and Zr II lines in the solar spectrum, we found that the non- 
LTE approach led to better agreement between the abundances derived from the lines of the 
two ionization stages than the LTE one. We could not constrain the scaling factor ku to the 
classical Drawin formalism due to the large spread in abundances derived from different lines. 
Nevertheless, we dwelt on the non-LTE case with allowance made for the hydrogen collisions in 
both permitted and forbidden transitions with kn = 0.1. The absolute zirconium abundance 
in the solar atmosphere (averaged over Zr I and Zr II lines) is \ogezr = 2.63±0.07. 
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Erratum: Formation of Zr I and II lines under non-LTE conditions of stellar 
atmospheres [Astron. Letters, 36, 664 (2010)] 

In the published paper, an outdated value of E ion (Zi I) = 6.84 eV was employed for an ioniza- 
tion energy of the neutral zirconium, instead of E ion (Zr I) = 6.63 eV measured by Hackett et al. 
(Hackett, P. A., Humphries M. R., Mitchell S. A., & Rayner D. M., 1986, J. Chem. Phys. 85, 
3194) and recommended by the NIST atomic data base (http:/ /physics.nist.gov/PhysRefData). 
The calculations with Ei on (Zr I) = 6.63 eV lead to lower number density of Zr I and, therefore, 
higher element abundance derived from the Zr I lines. For the Sun, such an abundance cor- 
rection amounts to 0.18 dex independent of the line-formation assumptions used. The change 
in E ion (Zr I) only weakly affects the departure coefficients of the important Zr I levels in the 
solar atmosphere, i.e., by no more than 1.5/3.2% inside logrsooo = —1/ — 2. This is because 
the ratio J v j ' B V {T) changes only a little when a wavelength of the ultraviolet ionizing photons 
increases by approximately 60 A. The change in Ei on (Zr I) does not affect the results for Zr II. 

The corrections to the published paper are as follows. 

1. Table 1, Table 3, and Fig. 7: the Zr I based abundances have to be increased by 0.18 dex. 

2. Abstract, Sect. 3.2, last paragraph, and Sect. 5, last paragraph. The revised abundance 
difference (loge(Zr I) - loge(Zr II)) is —0.10 dex at the LTE assumption and +0.10 up to 
+0.15 dex in the non-LTE calculations depending on the treatment of hydrogenic collisions. 
The calculations with E ion (Zr I) = 6.63 eV support our earlier conclusion that the Zr I/Zr II 
ionization equilibrium in the solar atmosphere cannot be established under the LTE conditions. 
The magnitude of the departures from LTE for Zr I can be overestimated due to missing high- 
excitation levels of Zr I in the model atom. 

The authors thank Nicolas Grevesse for picking our mistake up. 



